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Abstract 

The next generation of high energy e + e~ linear colliders is expected to operate 
a t y/s ^ 500 GeV. In this energy range the WW fusion channel dominates the 
Higgs boson production cross section e + e~ — ► uuh? /H°. We calculate the one-loop 
corrections to this process due to fermion and sfermion loops within the MSSM. We 
perform a detailed numerical analysis of the total cross section and the distributions 
of the rapidity, the transverse momentum and the production angle of the Higgs 
boson. The fermion-sfermion correction is substantial being of the order of —10% 
and is dominated by the fermion loops. In addition, we explore the possibility of 
polarized e + /e~ beams. In the so-called "intense coupling" scenario the production 
of the heavy Higgs boson H° is also discussed. 



1 Introduction 



The Standard Model (SM) of fundamental particles has been tested with an impressive 
precision by a large number of experiments. The resulting body of data is consistent 
with the matter content and gauge interactions of the SM and a Higgs boson h° of mass 
m h o < 204 GeV [1]. The four experiments at LEP delivered a lower bound for the SM 
Higgs boson mass, m h o > 114 GeV [2]. If a fundamental Higgs boson exists, it would fit 
very naturally into supersymmetric (SUSY) extensions of the SM, in particular into the 
Minimal Superymmetric Standard Model (MSSM). The latter requires the existence of 
two isodoublets of scalar Higgs fields, implying three neutral Higgs bosons, two CP-even 
bosons h°, H°, one CP-odd A , and two charged Higgs bosons H ± . The lightest Higgs 
particle h° could exhibit properties similar to those of the SM Higgs boson. Its mass 
is predicted to be less than 135 GeV [3], taking into account radiative corrections. The 
present experimental bound from LEP are m h o > 88.3 GeV and m A o > 88.4 GeV at 95% 
CL [2]. 

The next step in the search for the Higgs boson will take place at the Tevatron [4, 
5] in pp collisions at 2 TeV. The gluon-gluon fusion process is the dominant neutral 
Higgs production mechanism, but suffers from the overwhelming QCD background of bb 
production. The most promising Higgs discovery mechanism for m^o < 130 GeV is most 
likely the Higgsstrahlung qq — > W — > Wh°. The WW fusion process WW — > h°, i.e. 
pp — > q WW q —> qh° q, plays a less important role. 

At LHC, in pp collisions at 14 TeV, the gluon-gluon fusion mechanism provides the 
dominant contribution to Higgs boson production [6]. The next important Higgs pro- 
duction channel is the vector boson fusion VV — > h°/H°. In particular, it provides an 
additional event signature due to the two energetic forward jets. It has been argued that 
the channels WW — > h°/H° -> rf and WW can serve as suitable search channels at 
LHC, even for a Higgs boson mass of rrih ~ 120 GeV [7]. Very recently, it has been shown 
that the WW fusion process qq — > qh° q with h° — > bb may be used to identify and study 
a light Higgs boson at the LHC due the two rapidity gaps in the final state [8] 

The next generation of high energy e + e~ linear colliders is expected to operate in the 
energy range of y/s = 300 - 1000 GeV (JLC, NLC, TESLA) [9, 10, 11]. The possibility 
of a multi-TeV linear collider with y/s ~ 3 TeV (CLIC) is also under study [12]. At these 
colliders high-precision analyses of the Higgs boson will be possible. In e + e~ collision, for 
energies > 200 GeV, the production of a single Higgs boson plus missing energy starts 
to be dominated by WW fusion [13, 14, 15], that is e + e~ -> v e v e WW -> u e u e h°/H°, 



1 



whereas the Higgsstrahlung process [16] e + e~ — > Z/i — > z/z//i° becomes less important. 
The rates for the fusion are generally one order of magnitude smaller than those of 
the WW channel. 

The process e + e~ — > u e u e WW —> v e v e h° / H° was calculated at tree-level in Refs. [13, 
14, 15]. The leading one-loop corrections to the WWh° vertex in the SM were also calcu- 
lated (see the review article [17] and the references therein). For this coupling also QCD 
corrections were included, the 0(a s G F mf) corrections in Ref. [18] and the 0(a 2 s G F m^) 
ones in Ref. [19]. 

In this paper, we have calculated the one- loop corrections to the WWh°/H° vertex in 
the MSSM due to fermion/sfermion loops. We have also included the corresponding wave- 
function corrections to the W and Higgs bosons. They are supposed to be the dominant 
corrections due to the Yukawa couplings involved. We have applied the corrections to the 
single Higgs boson production in e + e~ annihilation in the energy range y/s = 0.5 — 3 TeV, 
i.e. to e + e~ — > v e v e WW — > v e v e h°/H°. We have also included the Higgsstrahlung process 
e + e~ — > Zh° — > vish°/H° and the interference between those two mechanisms. Because 
the Higgsstrahlung process is much smaller in this range, we have neglected its radiative 
corrections [20]. 

In a previous paper [21] we have already given some results for the light Higgs bo- 
son production. This work represents a much more detailed study of single Higgs boson 
(h° and H°) production in e + e~ collisions, including radiative corrections. In particular, 
important kinematical distributions of the rapidity, transverse momentum and the pro- 
duction angle of the Higgs boson are given. Polarization of the incoming e + /e~ beams is 
also included. In addition, we consider the "intense coupling regime" [22], where all Higgs 
bosons of the MSSM are rather light, and for large tan (3 couple maximally to electroweak 
gauge bosons and strongly to the third generation fermions. The paper contains a brief 
discussion on the background, although it has not been our intention to perform a Monte 
Carlo study. 

The paper is organized as follows: In section 2 we give the formulae for the tree-level 
amplitude, then calculate the one-loop corrections, taking into account fermion/sfermion 
loops. We express the corrections to the vertex in terms of form factors. In section 3 we are 
discussing the calculation of the cross section for the Higgs boson production, especially 
including the one-loop correction. In section 4, we perform a detailed numerical analysis 
and discuss our results. Finally, in section 5 we present our conclusions. Appendix A 
exhibits explicitly the expressions for the form factors. Appendix B gives details of the 
calculation of the cross section and the various distributions. 
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Figure 1: The Feynman graphs for the process e + e — > vi/h°/H°. For the Higgsstrahlung 
contribution |.M£ ee | 2 one has to sum over all three neutrino types. 

2 Matrix elements and one-loop corrections 

We study the process 

e + (p 2 ) + e~(pi) -> # fc °(p) + z>(p 4 ) + i>(p 3 ) • (1) 

p, pi, J92, P3, £>4 are the corresponding four momenta and = {h°, H }. 

The contributing Feynman graphs are shown in Fig. 1. The amplitude of Eq. (1) con- 
sists of three parts: WW fusion at tree level .M^ 6 , its one-loop correction AA4 W and the 
Higgsstrahlung process Mz ee , e + e~ — > Zh° — > P^i/^, i.e. .M = -M^ ee + .M| ee + AA^w- 
As already mentioned in the introduction, we have neglected the radiative corrections 
to A^^ ee , as this amplitude is much smaller than .M^ 6 in the energy region considered. 
We will include polarization of the incoming electron and positron beams. V- and V + 
denote the polarization of the e~ and e + beams, with the convention V± = { — 1, 0, +1} for 
{left-polarized, unpolarized, right-polarized} e ± beams, respectively, (e. g., V- = —0.8 
means that 80% of the electrons are left-polarized and 20% unpolarized.) The mass of 
the electron is negligible. Therefore, all vector particles propagate only transversally. We 
introduce the polarisation factors 

V LR = (1 - V-) (1 + V+) , V RL = (1 + V.) (1 - V+) • (2) 

Vlr (Pre) gets maximal in e\ (e^e£) collisions. 

The squared matrix element in the one-loop approximation is given by 

\M\ 2 = \MT? + \MT? + 2 3? [MT(M^ C ) ] + AM W {(M^ e ) ] + (Mf c ) ] )] . (3) 
The first three terms correspond to the tree level part [15]: 

|ACT = ^LRglwwd 4 JT2 _ m 2 YT7U _ m 2 n 2 Vx-V^Vi-Vz , (4) 



3 



o 4 1 
|A/ftree|2 _ 2 j/ ^ 

' Z ' ^ c 4 w (g 2 - m\f {{ql - m\f + m| T|) 

(C£) 2 P1.P4P2.P3 + Pi?L (C£) VP3P 2 .p 4 ) , (5) 

2$t[MT{MT) ] ] = 2V LR Clg kZZ g kWW ^ rPl .p,p 2 . P3 x 



?1 -m| 



(9? - ™!) ((d - O 2 + m\ r|) (*? - ro^fe* - m 2 ,) 



(6) 



k = 1,2, where 1 (2) stands for /i° (#°). 

The WW Hi couplings are <?iw^ = fl^w sin(/3 — a), g 2 ww = g m w cos(f3 — a) , gizz = 

mz sin(/3 — a), <?2zz = -^m z cos(j3 — a), f3 = arctan(t> 2 /t>i), a is the h°-H° mixing 
angle, k\ = p 3 —pi, k 2 = p 2 -p4, <?i = Pi +P2, q2 = P3+P4, is the total Z-boson width, 
C£ = sin 2 6 W — 1/2, C e R = sin 2 9 W: with 9 W being the Weinberg angle, and c w = cos 

Notice that the fusion is enhanced if the electron has left and positron right 

polarization. For instance, with V- = —0.85, V+ = 0.6, one has Vlr = 2.96 and Vrl = 
0.06. We are of course interested in the case where the fusion process dominates over the 
Higgsstrahlung process to get a large Higgs production rate. In this case the polarized 
cross section is just given by the unpolarized one times Vlr- 

Now we turn on discussing the calculation of the one-loop correction due to the fermion 
and sfermion loops. One expects them to be the most important corrections due to the 
Yukawa couplings involved. The renormalization of the five-point function simplifies to the 
renormalization of the WWH® vertex with off-shell vector bosons, where the renormaliza- 
tion of the other two vertices in the process (e. g. the e~u e W + coupling) is absorbed. The 
contributions of the first and second families of (s)fermions are numerically negligible due 
to the smallness of their Yukawa couplings. Therefore, we will consider the contribution 
arising from the third family of (s)fermions. 

At the one-loop level the Lagrangian for the WWH® coupling can be written as 

C = (g kW w gT + iV c (Ag kWW y») H° k W+ W~ . (7) 

The colour factor N c is 3 for (s) quarks and 1 for (s)leptons. Actually, for the calculation 
of the one-loop corrected WWH® vertex one has to compute the vertex and the wave- 
function corrections due to the graphs of Fig. 2, as well as the coupling correction Sg ( k c ^ vw , 

(Ag kWW y» = {5gtl w y v + (Sg^ + 5g ( k % w ) gT ■ (8) 

The vertex correction can be expressed in terms of all possible form factors, 

(5g k w w y U = + F lX KK + F 22 e 2 k v 2 + F 12 k^ 2 + F 2l e 2 k\ + 1 F^ s k lp k 25 , (9) 
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Figure 2: The Feynman graphs that contribute to the vertex (a)-(c) and the wave-function 
corrections (d)-(i). / (/') denotes the up (down) type fermion. 



k\ y 2 denote the four-momenta of the off-shell IV-bosons. At tree-level only the structure 
with g^ v is present, and therefore all form factors but F 00 have to be ultra violet (UV) 
finite without being renormalized. 
The wave-function correction is 



So (w) 

°ykWW 



9kww (| {8Zn)kk + SZw) + \ giww (8Zh, 



Ik 



(10) 



I ytz k, and 5Zh and 5Zw are the symmetrized Higgs boson and the TV-boson wave-function 
corrections calculated from the graphs (g)-(i) and (d)-(f) of the Fig. 2, respectively. 
In the case of the off-shell IV-bosons coupling to ei/ e , SZ\y has the form 



, 7 _ 5m 2 w - m^ w {kl) 5m 2 w - m^ w {kl) 5g 

O/jW — 5 1 g r £ 

w 9 



k\ — mw k 2 — m 2 



'111 
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The coupling correction is 

( c) (5g 5m w \ k sin 2(3 5 tan (3 

5 9*ww ~ [j + — ) 9*ww ~ (-1) —J^p9iww, (12) 

/ 7^ fc. The expressions on the right-hand sides of the Eqs. (10)— (12) can be found 
in Ref. [23]. Especially, we fixed the counter term 5 tan (3 by the on-shell condition 
SlU^m^) = 0, where IIaz^a) is the renormalized self-energy for the mixing of the 
pseudo-scalar Higgs boson A and Z-boson. By adding the vertex correction, Eq. (9), the 
wave-function, Eq. (10), and coupling correction, Eq. (12), we get the renormalized and 
therefore UV finite one-loop correction 

(Ag kWW y» = F°V" + F ll k±k\ + F 22 k%k v 2 + F 12 k^k u 2 + F 21 ^ + 1 F e e^ ps k lp k 25 , (13) 

which has exactly the same form as Eq. (9) but the form factor F 00 is substituted by the 
renormalized and hence UV finite one, 

pOO _ F 00 , r W , r (c) ( u) 

r — r + °ykww + °ykww ■ y 1 ^) 

Having calculated the form factors of Eq. (13), one can proceed to the calculation of 
the one-loop corrected cross section. The remaining parts of Eq. (3) due to the one-loop 
corrections are 



23? 



23? 



AM w (M%ry] = V LR g kWW g A (2 F°° Pl ■ p 4 p 2 ■ p 3 + F 21 s) 

(15) 

AM w (Mf c ) ] ] = V LR C e L g kZZ ^ 2 -(2F 00 Pl -p 4 p 2 -p 3 + F 21 s) 



X 



4 1 II <»> 



M - m\) M - m|) 2 + ml r|) 11 kf - m\, 



where 



S = (Pi ■ P4+ P2 ■ P3) (Pi ■ P2P3 ■ P4+ Pi ■ P4P2 ■ P3 ~ Pi ■ P3P2 ■ P4) 

- 2(p 1 -p 2 +p 3 -p 4 )p 1 -p A p 2 -p 3 . (17) 

As fci j2 are spacelike, F 00 and F 21 have no absorptive parts and are therefore real. The 
term with F € does not contribute to the cross section. The explicit forms of F 00 and F 21 
are given in Appendix A. 
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Figure 3: The integration variables in the rest frame of the two final fermions. 

3 Calculation of the cross section 

In order to calculate the cross section, including the radiative corrections from fermion and 
sfermion loops, one has to choose an appropriate reference frame. A detailed discussion 
on this is given in Appendix B. The momenta of the particles participating in the process 
are defined in Eq. (1). 

The calculation of the total cross section for the Higgs boson production e + e~ — > vvH® 
is performed in two steps. First, we calculate the differential cross section 



where E p is the energy of the produced Higgs boson. For this calculation it is convenient 
to work in the rest frame of the two final fermions, where one has ps + p^ = 0, see Fig. 3. 
In this frame the differential cross section can be evaluated using 



by integrating the total amplitude \ Ai\ 2 over the angles 6, <f), as they are defined in Fig. 3. 
For the tree-level case, it has been shown [14, 15] that these integrations can be performed 
analytically. For example, the results for the fusion process, the Higgsstrahlung, and their 
interference term from Eqs. (4)-(6) can be found in Eqs. (5)-(8) of Ref. [15]. One major 
complication of the inclusion of the one-loop corrections of Eqs. (15) and (16) is that it is 
not possible anymore to calculate these integrals analytically. This is due to the fact that 
the form factors F 00 and F 21 are functions of the momentum transfer k\ 2 - Therefore, for 
the one-loop corrected cross section we are bound to use numerical methods for this task. 




(18) 




(19) 
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Figure 4: The momenta in the the rest frame of the two initial fermions. 

The second step consists of the integration of the differential cross section in order to 
get the total for the Higgs boson production. To do this, we are working in the rest frame 
of the initial fermions, where p\ + f>2 = 0, see Fig. 4. In this reference frame we obtain 
the total cross section 



where 6 P denotes the angle of the produced Higgs boson with respect to the beam direction. 
Alternatively, one can use the rapidity y and the transverse momentum pt of the Higgs 
boson and calculate the cross section as 



where the integrand and the integration limits are given in Eqs. (45) and (46), respectively. 
These integrations, even in the tree-level case, are carried out numerically. The advan- 
tage of using the y, pt variables is the faster numerical convergence of the integration 
routines, due to the strong forward-backward peaking of the differential cross section for 
large y/s. For the calculation of the one-loop corrected cross section, one has to perform 
four numerical integrations successively. For this purpose, we use appropriate numerical 
integration routines found in the NAG library. In addition, we have checked that for the 
tree-level case our completely numerical calculation agrees with high accuracy with the 
semi-analytical results of Ref. [15]. 




(20) 
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4 Discussion and results 



In our numerical analysis, we have taken into account the contribution arising from the 
third family of fermion/sfermion loops. This contribution turns out to be the dominant 
one, in comparison with the first two families corrections, due to the large values of the 
Yukawa couplings h t and hf,. The impact of the running of the electroweak couplings g 
and g' is not negligible, especially for large y/s we are discussing here. Therefore it has 
been taken into account. For the calculation of the SUSY Higgs boson spectrum and the 
Higgs mixing angle a, a computer program based on Ref. [24] has been used. We note 
that for values of tan j3 > 5 and itla large the SUSY WWh° coupling mimics the SM one, 
while the WWH° is very small. This is due to fact that for these values of tan (3 we have 
sin(/9 — a) ~ 1 and cos(/3 — a) ~ 0. 

In the so-called "intense coupling regime" [22], where the neutral Higgs bosons are 
almost degenerate and light, tuh ~ rrih ~ m>A ~ 100 GeV, there is the possibility that 
the WWh° coupling is suppressed, while the WWH° one is not suppressed. For this case 
it is worth to explore the possibility of the heavy Higgs production. 

For simplicity, for all plots we have used A t — Aj, — A T — A, {mfj^m^^mi^m^} = 
{y|j, ^j, 1, 1} rag and M 1 = | tan 2 6 W M 2 ~ 0.5 M 2 . The choice of a common trilinear 
coupling and the correlation between the soft sfermion and gauginos masses are inspired 
by unification. 

Concerning the polarization, from the Eqs. (4)-(6),(15) and (16), we can see that 
basically the polarized cross section for the process e + e~ — > vvhP / H° is the unpolarized 
multiplied with the factor Vlr- Although the second term of the Higgsstrahlung contri- 
bution in Eq. (5) is multiplied with Vrl, considering that for y/s > 500 GeV the total 
tree-level cross section is dominated by the WW fusion channel, this term is not impor- 
tant for the Higgs production at future linear colliders. Therefore, if one wants to enhance 
the cross section the appropriate mode would be e^, where we have a po \ ~ Vlr Ounpoi, 
while Vlr being as high as 3 to 4. 

We will start discussing the light Higgs boson production e + e~ — > vvh , in the MSSM, 
and the impact of the fermion/sfermion corrections calculated in section 2. In Fig. 5 (left) 
we have plotted the total tree-level cross section (dashed line) and the one-loop corrected 
one (solid line) for up to 3 TeV. The tree-level cross section includes the contribution 
from the WW fusion channel, the Higgsstrahlung and their interference. The correction 
stemming from the fermion/sfermion loops is always negative and substantial being of the 
order of —10%. Focusing for y^i up to 1 TeV (right) in a more detailed figure, the various 
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Figure 5: The tree- level cross section <7o (dashed line) and the one- loop corrected one a 
(solid line) for y/s up to 3 TeV (left). The SUSY parameters are chosen as tan/? = 40, 
fi = -300 GeV, A = -100 GeV, m Q = 300 GeV, M A = 500 GeV and M 2 = 400 GeV. 
Focusing for y/s up to 1 TeV, the various contributions to the tree level cross section 
are presented (right). The dotted-dashed line represents the tree-level cross section o^ w , 
the dotted line the a^ w + a^~ str . The dashed line includes also the interference term 
0-mterf. anc j re p resen ts the total tree-level cross section. The solid line includes the one-loop 
correction. 

contributions are presented. The dotted-dashed line represents the contribution from 
the WW channel at tree- level alone, whereas the dotted line includes the Higgsstrahlung 
contribution as well. The dashed line includes in addition the interference between the 
WW channel and Higgsstrahlung. The size of this interference term is extremely small, 
hence the difference between the dotted and dashed lines is rather tiny. It can also be seen 
that for y/s > 500 GeV the WW fusion contribution dominates the total cross section 
for the Higgs production. Actually, for y/s > 800 GeV the total tree-level cross section is 
due to WW fusion. The solid line represents again the one-loop corrected cross section. 
In Fig. 5 we have taken: tan/? = 40, fi = -300 GeV, A = -100 GeV, = 300 GeV, 
Ma = 500 GeV, and M 2 = 400 GeV. Choosing different sets of parameters, the basic 
characteristics of these plots remain unchanged. Actually, the soft gaugino masses Mi j2 
affect only the Higgs boson masses and couplings through radiative corrections. 

In Fig. 6 the relative correction Aa/ao is presented as a function of y/s for two different 
sets of parameters. The solid line corresponds to the set used in the Fig. 5, whereas for 
the dashed line we have taken tan/3 = 10, /i = —100 GeV and A = —500 GeV, keeping 
the rest of them unchanged. This figure shows that the size of the one-loop correction to 
the Higgs production cross section is practically constant for y/s > 500 GeV and weighs 
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Figure 6: The relative correction Aa/ao as a function of \fs (Aa = o — er ), where ctq 
is the tree-level and a the one-loop corrected cross section. The solid and dashed lines 
correspond to two different choices of the SUSY parameters, as described in the text. 

about —15%, almost for any choice of the SUSY parameters. This is a consequence 
of the dominance of the fermion-loop contribution over the one-loop corrections, and 
therefore the total correction is not very sensitive to the choice of the SUSY parameters. 
This behaviour of the one-loop correction will be discussed further after presenting the 
influence of the corrections on the various distributions. 

In Fig. 7 we present the distributions and d ^ e as a function of the rapidity 

y, the transverse momentum and cos# p , respectively. We have fixed y/s = 1 TeV. 
The dashed lines represent the tree-level case, while the solid lines the one-loop corrected 
one. The SUSY parameters are as in Fig. 5. Paying attention to the figures of the ^ 
and rfc ^g distributions, we see that the tree-level distributions are completely symmet- 
ric. We have checked numerically that the one-loop corresponding distributions are also 
symmetric up to differences of O(10~ 2 ) fb. For example, at a collider like TESLA with an 
integrated luminosity of 500 fb -1 such a small asymmetry yields only few events, making 
these measurements very difficult. The reason for such tiny asymmetries is the smallness 
of the form factor F 21 , which contributes to the one-loop corrections in Eqs. (15) and 
(16). The dominant one-loop correction results from the form factor F 00 , which having 
the same structure as the tree-level WWh° coupling, that is a correction to the g^ u term 
in Eq. (7), does not contribute to the asymmetry. 

This fact, in conjunction with the behaviour of the correction illustrated in Fig. 6, 
suggests that a handy approximation of these fermion/sfermion loop corrections might 
be possible [25]. The major complication in calculating the corrections from Eqs. (15) 
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Figure 7: The distributions and d ^ g as a function of the rapidity y, the transverse 

momentum p T and cos^p, respectively. Here y/s = 1 TeV. The dashed line (solid line) 
represents the tree-level (one-loop corrected) distribution. The SUSY parameters are as 
in Fig. 5. 

and (16) is the dependence of the form factors F 00 , F 21 on the momentum transfer of the 
fused H^-bosons k\ 2 . On the other hand, the dominant contribution from the integration 
of these terms over the phase-space arises for small values of kf 2 - Therefore, the essence 
of such an approximation will be to keep the dominant Yukawa terms from the form factor 
F 00 for k\ 2 ~ 0, which will be just a factor correction to the tree-level coupling WWh°. 
In the literature there is an effective approximation, where one only corrects the WWh° 
coupling [18]. Although the sign of this approximation is correct, it does not however 
account for the whole effect. 

As it has been discussed earlier, in the bulk of the SUSY parameter space the WWH° 
coupling is rather small, resulting in a small production cross section for the heavy CP- 
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Figure 8: The various cross sections for the heavy Higgs boson H° production in the 
"intense coupling regime". The meaning of the curves in the plots is as in Fig. 5. The 
SUSY parameters are chosen as tan/3 = 30, \i = M 2 = 350 GeV, A = 1000 GeV, 
rriQ = 1000 GeV, M A = 110 GeV. 

even Higgs boson H°. The situation can be reversed in the so-called "intense coupling 
regime". There the WWH° coupling can be significant, while the WWh° coupling be- 
comes smaller. For this case we show Fig. 8. In order to approach this case the SUSY 
parameters have been chosen as tan/3 = 30, fi = M 2 = 350 GeV, A = 1000 GeV, 
vhq = 1000 GeV, M A = 110 GeV. For such a choice h°, H° and A are almost degenerate 
and light with masses from 110 GeV to 120 GeV. The meaning of the various curves in 
the plots here is as the corresponding in the Fig. 5. Comparing Fig. 8 with Fig. 5 we see 
that the cross section for the H° production is smaller than the h° production. Yet, it is 
possible to tune the SUSY parameters in such a way to obtain values for the heavy Higgs 
production as large as for the light one. In any case, it seems that going to the "intense 
coupling regime" the task of discriminating between the two Higgs bosons becomes not 
trivial. 

Fig. 9 exhibits the percentage of the sfermion loops to the total one-loop correction as a 
function of tan/3 (left) and \x (right), for two different values of /x and tan/3, respectively, 
as shown in the figure. In the left (right) figure we have chosen A = 100 GeV (A = 
400 GeV). The rest of the SUSY parameters are: = 300 GeV, M A = 500 GeV, 
and M2 = 400 GeV. Here y/s has been fixed to 1 TeV. The grey area in the right 
figure is excluded due to the chargino mass bound. We see that the maximum value of 
order 10% can be achieved for large values of \x and tan/3. There, due to the significant 
mixing in the stop and sbottom sector, the contribution of stops and sbottoms in the 
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Figure 9: The percentage of the sfermions to the total one-loop correction as a function of 
tan/3 (left) and \i (right). The rest of the SUSY parameters have been fixed as described 
in the text. Here \fs = 1 TeV. The grey area in the right figure is excluded due to the 
chargino mass bound. 

loops is enhanced. For these values of SUSY parameters the sfermion masses approach 
their experimental lower bounds. But even there the dominant correction, at least 90% 
of the total correction, is due to the fermion loops. 

Finally, let us discuss the background to single Higgs boson production, Eq. (1). To 
the background several processes contribute. Single Z-boson production from WW fusion 
has a large cross section (~ 200 fb at y/s = 500 GeV) and a similar topology as the signal 
process but the invariant mass of the two jets would peak at mz- Tagging of b would 
improve the signal as the branching ratio of Z — > bb is only 15%. The angular distribution 
of the two jets would also be different due to the spin 1 of the Z boson. Double Z-boson 
production, e + e _ — > ZZ, where one Z decays into vv and the other into two jets (a ~ 
500 fb at y/s = 500 GeV) is another background. It can be reduced by cutting out 
the forward direction and measuring the invariant mass of the two jets. An important 
background is due to the process e + e~ — > e + W~v e through jW fusion (a ~ 4.5 pb), with 
a very low px electron being lost in the beam pipe. Again the invariant mass of the two 
jets from W would give a peak at mw, and above all 6-tagging would strongly reduce this 
background [26]. Another source for the background is due to e + e~ — > e + e~bb (via 77 
fusion), where the e + and e~ are emitted in the very forward direction thereby escaping 
detection. However, the significance of this background can be substantially reduced by 
making a cut in |cos#| of the outgoing bb pair, to eliminate the part, where the bb is 
emitted near the beam direction [27]. 
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5 Conclusions 

In this paper we have calculated the fermion/sfermion loop corrections to single Higgs 
boson production e + e~ — > VvhP / H° in the context of the MSSM. They are supposed to be 
the dominant radiative corrections due to the size of the Yukawa couplings. At the next 
generation of high energy e + e~ linear colliders, where \/s > 500 GeV , the WW fusion 
channel dominates the Higgs boson production cross section. We have also included the 
Higgsstrahlung process and its interference with the WW fusion. The one-loop correction 
to the cross section is negative and of the order of —10%, and is rather independent of y/s 
for y/s > 500 GeV. It is dominated by the fermion loops, usually being larger than 90% 
of the total correction. For the case of maximal mixing in the sfermion mass matrices, 
the contribution of the sfermion loops is enhanced, but nevertheless weighs less than 10% 
of the total one-loop correction. The possibility of having polarized e + /e~ beams is also 
explored. The study of the kinimatical distributions of the rapidity and the production 
angle of the Higgs boson shows that these loop corrections do not alter the symmetry of 
the tree-level distributions. In the bulk of the parameter space of the MSSM the WWH° 
coupling is suppressed, making the heavy Higgs boson production very difficult. Yet, 
going to the "intense coupling regime" there is a possibility of obtaining sizable values for 
this coupling. We have studied the heavy Higgs boson production, including the one-loop 
fermion/sfermion corrections, in this case. 

Note added in proof 

After submitting our paper, it was claimed in the paper by T. Hahn, S. Heinemeyer, 
and G. Weiglein, hep-ph/02 11204, that there is a discrepancy between their and our 
results. The difference, however, was resolved in a recent paper by A. Denner et al., 
hep-ph/0301189. The difference in the size of the radiative corrections is due to the use 
of a different charge renormalization scheme and different input parameters. With the 
same input mass parameters we get the same result for the corrected cross section. 
We thank A. Denner and T. Hahn et al. for comparing results. Especially, we want to 
thank A. Denner for correspondence and the strong effort in clarifying the situation. 
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Appendices 



A Form factors F 00 and F 21 

For convenience we are presenting the formulae for the (s)top-(s)bottom doublet, but 
actually they are valid for any (s)fermion doublet. The form factors F 00 and F 21 are 
given by 

= J^ 2 ( A Z + A 2 + A % + A %), (22) 

Note that A 21 - = 0. A™, A 2 ] v A™, and A 2 tb \ correspond to Fig. 2a, Af- and Af- to Fig. 2b, 
and Af- to Fig. 2c. 



3 9 



AZ = - J {cos a, sin a} ((2 m 2 t +2m 2 b -k 2 - k 2 2 ) C - 8 C 00 

+5 (^, m 2, m fe 2 ) + B {k 2 , m 2 , m 2 ) + 2 B (m 2 H o,m 2 , m 2 )) , (24) 

3 2 

A 2 « = o 9 ^ A cosa,sma}(c + C 1 + A(C 2 + C 12 + C 22 )), (25) 
2 raw smp V / 

with C.. = C„(k 2 , m 2 , k 2 , m 2 , m 2 , m 2 ), and {cosa, sina} corresponds to k — {1, 2}. 

k 

3 2 

A Z = -^^-{-sma,cosa}((2m 2 t +2m 2 b ~k 2 1 -k 2 2 )C -8C 00 
4mw cosp V 

+B (kl m 2 , m 2 b ) + B (k 2 2 , m 2 , m 2 b ) + 2 B (m 2 H o,m 2 b , m 2 )) , (26) 

3 2 

= o g ^ o {-sina,cosa}(Co + Ci+4(C 2 + Ci 2 + C 22 )), (27) 
2miy cosp V / 

with C.. = C ..(k 2 ,m 2 Q ,k\,m 2 ,ml,m 2 ), and {— sina, cosa} corresponds to k — {1, 2}. 

h 

Aab o 2 \ r / /^(t D* TDb pb /nrl , /ob pb pb pt p< /nr2 \ fOCA 

^3? - ~ Z 9 2-^ \yi3k K il K jl K ll K ll U ba + ^ijk^il^jl^ll^ll ^ba J > l/°J 

i,j,l=l,2 

with the couple of indices (a&) = (00) or (21), k = 1,2 denotes /i° and if , C 1 = 
C (k 2 ,m 2 rr ,k 2 ,m 2 ,m 2 ,m 2 ) and C 2 = C (k 2 ,m 2 ,k 2 ,m 2 ,m 2 ,m 2 ). The definition of 
the rotation matrices it!* and i? b , and the coupling matrices Gr*- fc and G^ fe are given in 
Ref. [23]. 

2 



»,J'=1,2 
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B Calculation of the cross section 



In general the differential cross section for the process of Eq. (1), when the colliding 
fermions /i, / 2 are massless, is given by 



\M 



2 



da = -J !— d<S> 3 , (30) 

4pi -p 2 



where the 3-body phase-space is 



Jib (O ^4.4/ . \ d3 P3 d3 P4 d3 P /oi\ 

d*3 = (27T) ^(p.+^-^-p.-p)^^^-^^^. (31) 

The differential cross section can be cast into the form 

d 3 a f \M\ 2 4 d 3 p 3 d 3 p 4 , qo ^ 

Ep ^ = J T67(2^ 6{pi + P2 - p3 - p4 - p) -W -eT> (32) 

where s = (pi + p 2 ) 2 = 2p 1 ■ p 2 . 

In order to calculate the differential cross section of Eq. (32) we are following the 
procedure of Ref. [14] and we choose to work in the rest frame of the two final fermions 
defined by p 3 + p 4 = 0, see Fig 3. In this frame p — pi + P2, which means that the vectors 
p, pi and p 2 lie in the same plane, the (x, z) plane in our case. 

In this frame one gets 

= J T ^^t(Ei + E 2 -2E 3 -E p )dE 3 dcosed<j ) 

r 2n \m\ 2 

dcos9 / d<j> ' ' . (33) 
Jo s4tt 5 



The products Pi-Pj, i,j = 1, 2, 3, 4, which are involved in |.M| 2 , can be expressed in terms 
of the angles 9, 0, x an d the products p ■ pi, p ■ p 2 : 



Pi 


■P2 = 


s 

2 ' 




Pi 


■P3 = 


\{s-2p- 


Pi) (1 - cos6») , 


Pi 


•p 4 = 


\{s-2p- 


Pl) (1 + COS0), 


P2 


■P3 = 


\{s-2p- 


p 2 ) (1 — cos x cos — sin x sin cos 0) , 


P2 


■P4 = 


\{s-2p- 


p 2 ) (1 + cos x cos + sin x sin 9 cos 0) , 


P3 


■P4 = 


\(s + m 2 H 


-2p-p! -2p-p 2 ), 




kl = 


-2 Pi "PS, 






k 2 = 


-2p 2 -p 4 . 





(34) 
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It is important to note that cos x can be expressed in terms of Lorentz invariant quantities 

2s(s + m 2 H -2p-px -2p-p 2 ) 

cosx=l r . 35 

(s-2p-p 1 ) (s-2p-p 2 ) 

m H denotes the mass of the Higgs boson h°/H°. 

Using Eqs. (34), (35) and integrating over the angles <fi and 9 in Eq. (33) one obtains 
the differential cross section E p ^- as a function of the Lorentz invariants: cosx, P'Pi and 
p ■ p 2 . In order to calculate the total cross section we go to the rest frame of the initial 
fermions. There we have Pi + p 2 = 0, see Fig. 4, and we have chosen the beam direction 
as the z-axis. The three momenta of the final state particles p, p% and p^ span a plane, 
and 9 P is the angle between the pi and p. In this reference frame one finds 

p-p l = — (E p - \p\ cos9 p ) , p-p 2 = —(E p +\p\ cos 6 P ) , 
\p\ = y/E* - m\ . (36) 
The total cross section is given by 



a = ^ dfJjcoeO p J o d\p\\p? Yp [E p —^ 

= 2tt J dcos6 p J dE p \p\[E p —y (37) 



The maximum value of \p\ 
is for cos 5 = 1, see Fig. 4, and this gives 



p r x = s —^- (38) 



as integration limit for Eq. (37). 

It is more convenient to use in Eq. (37) the rapidity y and the transverse momentum 
Pt of the produced Higgs boson, instead of the integration variables E p and cos6 p . The 
transverse and longitudinal momenta of the Higgs boson are 

p T = \p\ sin 9 P , p L = \p\ cos d p . (40) 

Defining the rapidity as 
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one finds 



p L = m T sinh y , p ■ p 1 = ^ m T e y , 

E p = m T coshy , p-p 2 = — m T e y , (42) 



where the transverse mass of the Higgs boson is defined as 



m T = yj m\ + p\ . (43) 

Changing the integration variables from cos 9 P ,E P to y,pr, the total cross section of 
Eq. (37) can be cast into the form 



rv+ r 
/ dy / 

Jy- JO 



(#> m " , „ / d 2 a 
\dydp T 



a = I dy I dp 2 T [ ) , (44) 



where 

/ d 2 a \ /.. d :i a 



\ dy dp\ ) \ v d 3 p 
For the integration limits we find 



n E P -77— • (45) 



y ± = ±ln^, iA)™=( * + m \ ) 2 -m* H . (46) 



van ' \ 2 \fs cosh y 

One can reverse the order of the integrations and to carry out first the the integration 
over the rapidity y. In order to do this, we have to inverse the function {px) max (y) in 
Eq. (46). By doing this and by studying the integrations limits in the plane {y,p\) one 
gets 

f(Pr) m " rv+ f d 2 a \ 

°-L **i *(**?)■ (47) 

For the new integration limits we find 

/ _ 2 \ 2 

' ~2 \max S ~ m H 



where 



(p 2 T ) max = ( J , V± = ±ln z + v^T , (48) 



S + ^ (49) 



2 a/s (Pt + mj) 
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